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Abstract 

A system of three electrolytes separated by two parallel planes is con- 
sidered. Each region is described by a dielectric constant and a Coulomb 
fluid in the Debye-Hiickel regime. In their book Dispersion Forces, Ma- 
hanty and Ninham have given the van der Waals free energy of this system. 
We rederive this free energy by a different method, using linear response 
theory and the electrostatic Maxwell stress tensor for obtaining the dis- 
persion force. 
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1 Introduction 



The modest aim of the present paper is to give another derivation for the van der 
Waals free energy (generating a dispersion force) in a macroscopic arrangement 
of electrolytes. The considered arrangement is as follows: Space is divided into 
three regions i — 1, 2, 3 by two parallel planes distant by I from each other; with 
the Cartesian coordinate system chosen such that the x axis is normal to the 
planes, we define region i — 1 by x < —(1/2), region i = 2 by —(1/2) < x < (1/2), 
and region i = 3 by x > (1/2). Each region may be filled by an electrolyte. The 
solvent is macroscopically described by a continuous dielectric medium with a 
given dielectric constant Cj and a dissolved salt microscopically modelled by 
a Coulomb fluid in the Debye-Hiickel regime with a given Debye wavenumber 
Ki. In their book P^, Mahanty and Ninham write, for this arrangement, a free 
energy, per unit area of one plane, of the form 

m =^fr w in[i ■ Ai2 Ai3 exp( - 2s2 ° ] ' a) 

where fcs is Boltzmann's constant, T the temperature, sf = k| + k 2 , and A y - = 
(e»Sj - 6jSj)/(eiSi + ejSj). If k 2 = 0, F(l) is long-ranged; if k 2 ^ 0, F(l) is 
"screened" into a short-range expression. The result has many applications 
in electrochemistry, biology, etc. 

In ref. £Q (see also pj]), the authors start with a quantum system. They work 
out a quantum elaborate formalism, in which the /-dependent part of the free 
energy is written as a sum on Matsubara frequencies = 2imk^T '/%, where 
Ti is Planck's constant (divided by 2n). Actually, eq.Q) is the £ n = term of 
that sum. In the high-temperature regime, the quantum free energy should go 
to the classical one, and only the = term JIJ should survive. Therefore, a 
derivation of l|T]l. using from the start equilibrium classical (i.e. non-quantum) 
statistical mechanics, should be possible. Such a derivation has been given 
by Netz in a field-theoretical formalism Here, we present another purely 
classical derivation, based on the use of linear response theory and the Maxwell 
stress tensor, generalizing a method already used in a special case |J; in this 
way, we obtain the dispersion force corresponding to QJ. 

In classical statistical mechanics, there is no coupling between charges and 
electromagnetic radiation (Bohr- van Leeuwen theorem) [SJEIEI- Therefore, we 
have to consider only charges interacting by the static Coulomb law. 

We model each electrolyte as a system of species a, = 1,2, ...,Mj of 
particles, each of them carrying a charge q ai . The fugacity of species on is given 
as z ai . If each electrolyte occupied the whole space, the number density of 
species ai would be n ai ■ The Debye wavenumber m in region i is defined by 

2 ^Ea, ™«i<4 (r> , 
K i = , I 2 ) 

where (3 = l/(k B T). 

The paper is organized as follows. In Section 2, we compute the electric 
potential created by an unit charge. In Section 3, we use the linear response 
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theory for computing the Maxwell electrostatic stress tensor. In Section 4, we 
compute the kinetic part of the stress tensor. Finally, in Section 5, we obtain 
the total stress, and, by integration, recover the Mahanty-Ninham free energy 
P. 

2 A Green function 

In the following, we shall need the electric potential at ri knowing that there is 
a unit charge at r2. In the Debye-Hiickel regime (which describes the small- j3 
case) this potential G(ri,r2) is the Green function obeying 

(A x -^)G(n,r 2 ) = -(47r/e i )*(r 1 - r 2 ), (3) 

where the index i denotes the region where is ri and j denotes the region where 
is Y2- At the lowest relevant order in (5, the densities can be replaced by the 
bulk ones, this is why the bulk Ki are used in ©. This equation is supplemented 
by the continuity conditions at the interfaces of G and e(dG/dxi), and by the 
condition that G should vanish when \xi \ — > oo (this last condition means that 
we assume that there is no bulk potential difference between regions 1 and 3, 
and these potentials have been chosen as zero without loss of generality). 

Eq. @ can be replaced by a simpler one for the transverse Fourier transform 
of G. Let r be the two-dimensional vector which describes the components 
of r normal to the x axis. Since G depends on and only through their 
difference, we can define their transverse Fourier transform by 

G(xi,x 2 ,k)= /d^-r^exphik-^-r^G^,^), (4) 

where k = |k| is the modulus of a two-dimensional wavenumber normal to the 
x axis. Eq.© becomes a simpler differential equation. We shall only need the 
case when X2 is in region 2. Then, this equation is 

d 2 

(7^2 - fc 2 - n 2 )G(x 1 ,x 2 ,k) = -(4n/e 2 )S(x 1 - x 2 ). (5) 
The solution of this equation which vanishes when \x\ \ —> 00 is of the form 
G = cxp(-s 2 |a;i — X2V) + A(x2) exp(-s 2 £i) 

£2^2 

+S(a; 2 )exp(s2a:i) when - (1/2) < xi < (1/2), 

G = C(x2) exp(sixi) when x\ < — (1/2), (6) 

G = D(x2) exp(— S3X1) when x\ > (1/2), 

where sf = k 2 + k 2 . The coefficients A,B,C,D are to be determined by the 
continuity conditions at the interfaces. One finds, when both x\ and X2 are in 
region 2, G = Gbuik + G\, where Gbuik(^i, x 2 , k) = (27r/e 2 S2) exp(-s 2 |a:i - a^l) 
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is the Fourier transform of the bulk Debye screened potential Gbuik( r i, ra) 
(exp(— K 2 \ri — r2|)/(e2|ri — r2 1) and Gi, the Fourier transform of G\, is 

- 2n 2A 12 A 32 exp(-2s 2 l)cosh[s2(x 1 - x 2 )} 

Gi{xi,x 2 ,k) = — 

e 2 s 2 1 - A12A32 exp(-2s 2 t) 

2tt exp(-s 2 /){Ai2 exp[-s 2 (xi + x 2 )] + A 32 exp[s 2 (xi + x 2 )}} 

e 2 s 2 1 - A12 A 32 exp(— 2s 2 l) 



(7) 



3 Maxwell stress tensor 

In statistical mechanics qG(ri,r 2 ) means the statistical average of the total 
electric potential at ri when an additional charge q is placed at r 2 . Therefore, in 
terms of the microscopic electric potential </>(r) created by the system (dielectrics 
plus electrolytes, NOT q), 

< 0(n) >,= <zG(n, r 2 ) - q (8) 

Fi - r 2 | 

where < . . . > g is a statistical average taken when g is present. The presence 
of q adds to the Hamiltonian the term H' = q(f>(r 2 ). In the classical regime, 
assumed to be valid for both the charged particles and the dielectrics, linear 
response theory says that, when q is infinitesimal 

<4>(T X )> q =-P<H'4>{Tx)>, (9) 

where < . . . > is a statistical average taken without H 1 in the statistical weight. 
In writing ©, we have taken into account that < 4>(yi) >— at the lowest 
relevant order in (3; indeed, if the rhs of eq.© is replaced by zero, its solution 
vanishing at infinity is identically zero. We shall use JHJ and for computing 
the /-dependent part of the Maxwell stress tensor. We define a /-dependent part 
by the condition that it vanishes when / — > 00. For obtaining that /-dependent 
part, it can be checked that it is sufficient to replace the rhs of iJSJ by its part 
gGi(r!, r 2 ). From now on, all quantities refer to their /-dependent part. (JSJ) and 
(0 give 

(3 < ^(ri)^(r 2 ) >= -Gifa.ra). (10) 
The xx component of the Maxwell stress tensor 1 in region 2 is [S] 

T xx (*i) = £ (g^- 2 - \ V i ' < <K r i)0(r 2 ) >r 2=ri , (11) 
or, in terms of the Fourier transform of {TUJ), 

W (ft^ - k2 ) G^ u x 2 ,k) x „ 

1 Since here the electric potential is the microscopic one due to the charged particles of 
the system, not the macroscopic one due to some external charges, the simple Minkowski 
expression of the stress tensor is appropriate 
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dk 



{k\ + 2fc 2 )2Ai 2 A 32 exp(— 2s 2 /) 



(27r) 2 s 2 [ 1 - A12 A 32 exp(-2s 2 /) 



k 2 exp(— s 2 Z) [A12 exp(-2s 2 x) + A 32 exp(2s 2 a;) 



1 - A12 A 32 exp(-2s 2 /) 



(12) 



4 Kinetic part of the stress tensor 

The xx component of the stress tensor has, in addition to its Maxwell part, a 
kinetic part, which is, in region 2, — &sTn 2 (a:) where n 2 (x) is the total number 
density of the electrolyte at x. It happens that the departure of n 2 (x) from its 
bulk value n 2 = X) Q2 n ai gives to the /-dependent part of the stress tensor a 
contribution of the same order than the Maxwell part, although the Maxwell 
part can be and has been computed, at the lowest order in (3, by using the bulk 
densities in 

For computing n 2 (x), a natural idea is to use the first equation of the YBG 
hierarchy which gives dn 2 (x) / dx; however, it is not easy to find the integration 
constant for getting n 2 {x). We resort to a slightly more complicated method 
relying on the diagrammatic renormalized Mayer expansion reviewed in ref. [5]. 
The Debye-Hiickel regime is described by writing for the excess free energy F exc 
only the diagrams 



ring diagrams , 



(13) 



where a black circle stands for the density of some species of particles at some 
position (to be summed and integrated on), a bond stands for — (3 times the 
Coulomb interaction between the particles at its ends, and the ring diagrams 
are rings made of n (n > 2) black circles and n bonds. Taking the functional 
derivative of (fT3"|l with respect to n a2 (r) gives 



In 



o • 



+ ...), (14) 



where all white circles refer to the same fixed particle of species a 2 at some point 
r in region 2 (actually the densities depend only on the x component of r). In 
the rhs of (|14|) . the first diagram is —f3q a2 < 4>{x) >. The other diagrams are 

the sum of chain diagrams minus the first one o o; — (3 times the screened 

potential q\ 2 G is known to be given by this sum of chain diagrams (taking into 
account the position-dependence of the densities in the calculation of G would 
give corrections of higher order in /3). Therefore, 



In 



(^) 



-/3q a2 < <j>(x) > 



lim 



G(r,r') 



1 



£2|r - r'| 



(15) 



For the bulk values, a similar equation is obtained, by replacing in (|15fl n a2 (x) by 
n a2 , < <j){x) > by < 0buik >, arid G by Gbuik- Combining these two equations, 
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we obtain 

n a2 (x) = n a2 exp 



-0q aa [< 4>{x) > - < 



Gi(r,r) 



(16) 



Linearizing the exponential in Q15[l. summing on «2, taking into account the 
bulk neutrality condition y] q a2 n a 2 = which supresses the term containing 
< 4>{x) > — < 0buik > 7 and using J7J, we obtain 



n 2 (x) -n 2 = — — 



dk 



(2^) 2 s 2 



2Ai 2 A 32 exp(-2s 2 Q 
1 - A12 A 32 exp(-2s 2 Z) 



cxp(— s 2 /)[Ai 2 exp(-2s 2 x) + A 32 exp(2s 2 Z)] 



1 - A i2 A 32 exp(-2s 2 Z) 
Since n 2 (x) — n 2 — > as I — > oo, the Z-dependent part of ni(x) is n 2 (x) — n 2 



(17) 



5 Final result 

Finally, by adding the Maxwell and kinetic contributions in region 2, one finds, 
for the Z-dependent part of the force per unit area acting on the region on the 
right of x, f = -T xx (x) + k B T[n 2 (x) - n 2 ]. Thus, 

k B T r dk 25 2 A 12 A 32 exp(-2 S2 Q 
/_ 2 J (2tt) 2 1- A 12 A 32 exp(-2 S2 0" 1 J 

This force is, as it should be, independent of x. This general expression has 
been rederived in special cases, sometimes by different methods, see for instance 
ref. |101 B|. Since / = —dF(l)/dl, integration of / with respect to I, with the 
integration constant determined by the condition that F vanishes as I — > oo, 
gives Q. 

In the formalism used by Netz [3], Q has been obtained in a "one-loop" 
approximation; a "two-loop" correction has been computed in a special case. 
An alternative approach for obtaining corrections would be to extend the cal- 
culation in the present paper by taking into account higher-order terms in the 
renormalized Mayer expansion. We do not know which approach is the simplest 
one. 
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